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Abstract 

The emission of energetic gluons in tt production in e + e~ annihilation can 
have important experimental consequences, in particular on top quark mass 
measurements. We present compact, analytical expressions for the gluon energy 
distribution and its average value at first order in QCD perturbation theory. Our 
results are valid for arbitrary masses, collision energies and production currents. 
We pay particular attention to top quark production near threshold, and show 
that in certain cases the soft gluon approximation is insufficient to describe the 
radiation spectrum. 
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1 Introduction 



The discovery of the top quark - one of the basic components of the Standard Model 

- is one of the most important goals for present and future experiments. Indirect 
evidence for the existence of top is very strong, see for example |IJ, and the chances 
that it will be detected at the Fermilab pp collider in the next few years are very high. 
Once top is discovered, the next challenge for experiments is to measure its parameters 

- in particular its mass M - as precisely as possible. 

An unambiguous interpretation of experimental data and the determination of the 
top quark parameters relies on a clear quantitative understanding of the details of the 
production process, including the effects of gluon bremsstrahlung at the production 
stage [0, [|. To quantify these effects, one can study the average amount of energy lost 
by the top quark - the simplest infrared-safe characteristic of gluon bremsstrahlung. 
Such a quantity is of potential practical importance in top mass determinations and, 
for energies above threshold, can be calculated quite straightforwardly without undue 
complications from non-perturbative phenomena or top- width effects [0, E3, M. The key 
point here is the large top width, T t fj. Since we know from experiment that the top 
mass is greater than about 120 GeV |J, the decay width is greater than the typical 
hadronic scale, r t > // ~ 1 fm" 1 , and so the top quark decays before non-perturbative 
effects like the formation of t- flavoured hadrons become important 0, [7J . 

The bremsstrahlung of gluons off the top quark can have an important impact 
("radiation damage") on various characterstics of the tt final state. Taking e + e~ — > tt 
as an example, it obviously affects the kinematical constraint E t = -Ebeam, over and 
above the well-known effects of initial radiation. | It also influences the form of the tt 
production vertex and measurements of the spectra of secondary leptons from t-decay 
0,0. 

Since one of the primary goals of a future e + e collider would be to obtain an 
accurate measurement of the top mass M, it is worth considering the impact of gluon 
radiation on such a measurement. There are two scenarios which have been considered 



12]. If the top is not too heavy, it can have sizeable kinetic energy and a mass 
measurement can be obtained from the invariant mass per hemisphere. In practice, 
the distributions in the heavy (Mjj) and light (Ml) hemisphere masses can be fitted to 
a Monte Carlo simulation with M as a free parameter. However, a gluon radiated from 
the tt system will distort this measurement. In particular the average difference in the 
hemisphere masses squared is proportional (at 0(a s )) to the average gluon energy, 

{M% - Ml) = yfi{E g ) . (1) 

The more energy is carried away by the gluon, the more asymmetric the hemisphere 



2 These initial state radiation effects can be incorporated in a standard way |pL [| and will not be 
addressed further here. 
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masses. In practice, as we shall see, (E g ) can be 0(10 GeV), which is much larger 
than the estimated statistical errors on M from such analyses. A Monte Carlo which 
does not take full account of gluon radiation at the it production stage could yield 
misleading results on the precision on M. 

Another method which is more suited to a heavier top quark is to reconstruct the 
mass directly from the three-jet final state [|12|. To avoid combinatorial and other 
backgrounds, the energy of the three-jet system is typically constrained to be close to 
the beam energy. But once again, the emission before the top decay of an energetic 
gluon can distort this constraint, since naively 

(E3 jet) — -^bcam = —\(E g ) . (2) 

Once again we see the importance of taking the effects of gluon radiation into account 
in reconstructing the top mass. 

In this paper we present results for gluon bremsstrahlung in e + e~ — > ti which are 
of practical interest for physics studies of future linear e + e~ colliders. We focus in 
particular on the gluon energy spectrum which, as we have seen, can have important 
implications for the top mass measurement. It is also straightforward, at least in 
principle, to extend our results to the more complicated case of hadronic ti production. 

Since centre-of-mass energies -y/s at which multiple gluon radiation should be taken 
into account |fj, |T3[ (log(i/i/M) 1) appear unrealistic for the foreseeable future, 
we need only consider a first-order perturbative analysis. This is calculationally quite 
straightforward, and in fact many parts of the calculation can be taken over from 
existing QED analyses [13], [L5| (see also fl6|j ). 



The remainder of the paper is organised as follows. We first set up the theoretical 
framework which allows us to calculate the gluon radiation distribution for arbitrary 
production currents. We derive compact, analytical expressions for the gluon energy 
spectrum and its average value for vector, axial vector, scalar and pseudoscalar pro- 
duciton currents. We then present some numerical calculations which illustrate the 
size of the effects for typical experimental parameters. Finally, we summarize our 
results and present our conclusions. 



2 Calculation of QCD bremsstrahlung distributions 

In this section we describe the calculation of the cross section for single primary gluon 
emission in the process e + e~ — > ti. Our calculation is valid for any collision enegy 
■s/s > 2M, and we take both photon and Z exchanges into account. 

We begin by defining some variables: \fq^ = y/s = IE is the total centre-of-mass 
energy and M is the mass of the top quark. The momenta of the initial and final 
state particles are labelled by t~(k\) + t r (k'i) — > t(pi) + i{j>2) + g(k) and we define 
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= Pi + p% + k^. The energy fractions shared by the t, t and gluon in the e + e 
centre-of-mass frame are 

2g-p» -19 

2 = — 5— , 2 = 2 - zi - z 2 ■ (3) 
Q 

The quark velocity in the tt centre-of-mass frame is 



P = P(z) = s]l - po = y 1 - ^ < u = V 1 -^ , (4) 

where 



7 = — < 7 > Po = ~i : • ( 5 ) 



M 2 1 4 7 
4 1 — 2 

The final-state radiation is limited by the maximum kinematically allowed energy 

< z < z max = v 2 . (6) 

The virtualities of the top quark and antiquark before gluon emission, Ki, can be 
written in terms of the energy fractions: 

K X = 2k-pi = q 2 (l - z 2 ) , 

k 2 = 2k-p 2 = q 2 (l-zi) . (7) 
It is also convenient to introduce the "angular" variable 

«i 1 - z 2 

P = — = • (8) 

zq z z 

In the calculations which follow, we neglect the electron mass m e . Since we are con- 
cerned with energy losses caused by gluon emission at the tt production stage, we may 
treat the t-quarks as stable objects. Only very soft gluons with energy E g £ T t E/M 
are affected by the top width, and these make a negligible contribution to the average 
energy loss.Q 

2.1 Lowest-order cross section 

The total non-radiative cross section for tt production in e + e~ collisions can be ex- 
pressed in terms of the vector and axial contributions (see for example |T~9f| ) : 



on = ° VV (s) + a AA (s) = Ru(s)a pt (s) , (9) 



3 A detailed discussion of the potential impact of the top quark instability on gluon emission can 
be found in Refs. g |, |, 0. 



3 



where 

M*) = -^r 1 (io) 

is the point cross section for e + e~ — ► [i + with the QED coupling evaluated at the 
JfcOO = + i^ A (s) = v (( V T V (s) + ( a t a (s)) (11) 

with 

(v = 1(3 -^ 2 ) = 1 + 2 7 , (12) 
(a = v 2 = l-4 1 . (13) 

The electroweak factors t c (C = A, V) are defined as 

v 4 „ sis -Ml) 



3 et (s- MlY + {sT z /M z y 

s 2 

.2 i „2\ 2.2 s 



with 



+3(Ve + K ( 8 - miy + ( 8 r z /M z )=» ' (14) 
rA = %2 + a ^' K2 F¥T(«' (15) 



V2G F M Z 
47ra(s) ' 
-| + 2sin 2 6' w , a e = — \ , 



v t = | - - sin 2 W , a t = \. (16) 



4 
3 

The functions t v '(s) and t a (s) are shown in Fig. 1. Except for energies close to the Z° 
peak, which are not relevant for ti production, the axial current piece is numerically 
small compared to the vector current piece. At very high energies the ratio of r A to t v 
tends to a constant value of approximately 0.26. Furthermore, close to ti threshold the 
axial current contribution is further suppressed by two powers of v, see Eqs. (|I~2|,|13|) . 
We should also note that near threshold, the ti cross section is strongly modified by 



Coulomb enhancement effects, see Refs. 20, Ell E2|. However the extreme threshold 



region is not of primary interest in the present study, since the perturbation radiation 
is strongly suppressed there. 

2.2 Dalitz plot distribution for the radiative process 

We consider in this section the calculation of the differential cross section for the 
emission of a primary gluon of momentum k in ti production. After integrating over 
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the relative orientation of the quark and gluon momenta with respect to the incoming 
leptons, the interference between the vector and axial pieces vanishes and the double 
differential cross section in the quark and antiquark energy fractions (the Dalitz plot 
distribution) can be written as 



d 2 W 1 d 2 a g R vv j d 2 W \ R AA f d 2 W 



with 



dz±dz 2 (Tu dz\dz2 Rtt 1 dz\dz2 J T/ Rtt 1 dz\dz 2 



d 2 W \ 1 d 2 a c g 



dz\dz2) c o cc dz\dz2 



(17) 



(18) 



The subscript C (= V,A) denotes the production current, see Eqs. (P,P1) . 

The derivation of the expressions for the radiative matrix element squared, summed 
over colours and spins and integrated over the azimuthal orientation of the produced 
particles with respect to the beam direction, drastically simplifies if one uses an ana- 
logue of the "invariant-integration technique" which was employed in Refs. [13], [15] for 
calculations of the QED radiation accompanying two-charged-particle production in 
e + e~ annihilation. 

The matrix element Mq describing the radiation off the top quark produced via 
the C-current exchange can be written as 



with 



M§> = 9s BcT a ^ {u(Pi)A°v(p 2 )} ■ , (19) 



^ = rJ (j .e A ) + ^ + &. (20) 

Here T a is a SU(3) colour matrix and e\ is a polarization vector for the gluon. The 
matrices (C = V, A) are 

r^ = 7 M , 1^ = 7^75. (21) 

The four- vector j M is the classical current induced by the acceleration of the top-quark 
colour charges, 

U = — - 2 —, (22) 
and J^ c describes the C-exchange contribution at the electron vertex, 

Jf = v^T^h) . (23) 
The normalization factors Br are chosen so that the non-radiative matrix element 



MP is 



M§> = B C 9 -^ {u( Pl )^v( P 2)} ■ r u ' c . (24) 
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Following Refs. [HI it is convenient to integrate out the angular variables 
which describe the orientation of the final state quarks and gluons with respect to the 
incoming leptons. We first introduce the Lorentz-invariant polarization tensor Uz u 



n 



c 

flU 



-£ Tr [fa + M)A°(M-p 2 )A° 

A 



(25) 



with = ■y°(A^) + 'j , which describes the radiative decay of a polarized spin-one 
heavy "object" C. Since we wish to integrate over the directions of the quark and 
gluon momenta, the resulting tensor structure can depend only on the four-vector 
n v = Qv/ which defines the reference frame in which the quark energies are fixed. 
It is convenient to write the polarization tensor Il£, in the general form 



n 



c 

flU 



c > 



(26) 



where 



n c = n 



it /j 



TT 9 



(27) 



Note that because of vector current conservation, the tensor UY V is transverse, i.e. 

K = • (28) 



For the axial current, the second (longitudinal) term in Eq. fl26|) vanishes in the limit 
m e — > 0. In what follows we will set m e = 0, and so only the transverse component 
He appears in the final expression for the double-differential distribution, Eq. (|18|). 
As a result, we obtain 



d 2 W 



where 



c 



■>c 



1 CpOis 
V 



71 



(29) 



i n 



c 



4fi( 0) 



with ng } the non-radiative polarization operator that one obtains from Eqs. (|26| , ^7D 
by substituting the Born vertices for A° in Eq. 



rt(0)c 



16q- 



Tr[(p 1 + M)rJ(M-p 2 )f 



With the numerical normalization chosen for Eq. (HE) one has 



411 



(0) 

c 



Co, 



S c = — U c 



(30) 
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where the Born factors Cc are given in Eqs. ([T^,|TB|). Explicit calculations giveQ 



with 



U v - 

n A = 



. 2 <r > 1 /«2 , Kl 
4 2 Vki k 2 



■2<f , 1 



(4 2 



«2 «1 
Kl «2 



;i + 2 7 ) + 4 7 



Kl _ K2 _ 2M 2 ) - M 2 ('^ + ^ 

«1 K 2 



(31) 
(32) 

(33) 



In terms of the dimensionless energy fractions defined in Eq. ([3]), we have 

z\ + z 2 - 1 - 27 7 7 



4 J 



„2 



^2 J 



^2 



(1-^)(1-^) 



- 2 . 



(34) 



Finally, the Dalitz plot distribution is given by substituting the expressions given in 
Eqs. (y,0,0,|l-|l into Eqs. (0©- 

Before presenting numerical results, we make some general comments concerning 
the global structure of the above results for the He- 



The first ("classical") term in Eqs. (31.32|) corresponds to long-distance radiation 
with polarization vector e\\ in the plane defined by pi,p 2 - It leads to a universal 
contribution to the radiative cross section. In contrast, the short-distance effects 
are represented by the second term, which induces a current-dependent, non- 
universal contribution. This describes equal production of the gluon states eji 
and e*j_ (polarization transverse to the pi,p 2 plane). 

Next, consider the case of soft gluon emission, when the radiative cross sec- 
tion can be expanded in powers of the gluon energy fraction z < 1. In accor- 
dance with the factorization properties of soft radiation [23|, the first term in 
Eqs. fl3"IU3"2p is proportional to the product of the Born term (c and the usual 
accompanying radiation factor j 2 , with the former evaluated at the total centre- 
of-mass energy q 2 of the process. It is interesting that at small z the second 
(non- universal) term in Eqs. ( |3~1"| , |3"2"|) appears to be 0(z 2 ) compared to the first 
term, not 0(z) as one might naively expect. The only correction linear in z is 
incorporated into the classical radiation factor j 2 . The physical origin of this can 
be understood p4|| by recalling the celebrated Low soft bremsstrahlung theorem 
25H, extended to the case of charged fermions bv Burnett and Kroll [EB1. The 



4 For completeness, we present in the Appendix analogous formulae for scalar (C 
doscalar (C = P) exchanges. 



S) and pseu- 
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application of these classical results is especially transparent for the case of the 
integrated annihilation cross sections, where the non-radiative process depends 
only on one energy variable (the centre-of-mass energy of the charged particles) 

(hi) In the extreme ultra- relativistic limit, 7 < 1, the Hy and Ha distributions 
become identical. 



(iv) The results (|29| - |32|) are not valid in the extreme threshold limit, even within the 
context of perturbation theory. In this region, the first-order radiation amplitude 
is strongly modified by the QCD Coulomb-like interaction between the t and i 
[pTH. This leads to a distortion of the radiation pattern and to the additional 



suppression of radiation [21, 



2.3 Gluon energy spectrum 

To calculate the inclusive gluon energy distribution it is convenient to rewrite the 
quantities in terms of the variables z and p (see Eqs 

1 



3(3-|32])). This gives 



S so ft + \z 2 



J soft 



±z 2 



1 

1 

u 



- 1 



p 

1 27 
p p 



+ (l<->2) 



2). 



(35) 
(36) 



Here we have introduced E so f t as the classical current contribution which accounts for 
the z° and z 1 terms, 

„ 1 — z — 27 7 



J soft 



p 



and 



i(l-/^))<p<|(l + /^)). 
Explicit integration then gives 



dz J 



z 



c 



a 



cc 



dz 



(3(z) C F a s 



(37) 
(38) 

(39) 



with the Born factors Cc given in Eqs. ([L2],[T3]) . The universal "soft" part of the gluon 
density reads 

£{z) = (1 + v 2 - 2z)C{z) - 2(1 - z) , (40) 

5 The Low-Kroll-Burnctt approach leads to simple results also for initial-state radiation in anni- 
hilation processes. For example, in e + e~ annihilation, the universal part of the integrated photon- 
emission cross section is given by the product of the Born term taken at the point (q — k) 2 and the 
radiation factor describing the emission off the initial electrons [|l4| |l5|, |24|] . This construction remains 
valid for the more general case of arbitrarily polarized, massive initial-state particles (see also pTj ) . 
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where 



C(z) 



1 ln .l+^) 



(3{z) 1 - (3(z) ' 

The "hard" contributions which depend on the the production channel are 

Zh(z;V) = C(z)-1, 



Zh(z;A) = C{z) 



3-v 2 



1 . 



Note that the result for j^j}^ coincides with that obtained in Ref. [[14]]. 
2.4 Energy integrals 

We turn now to the quantity which is the main concern of the present study 
average energy lost by the i-quark by gluon emission at the production stage. 
The average gluon energy fraction is 



(z)c 



dz z 



dW\ 

. dz J 



c 



and the mean energy fraction lost by the t-quark (antiquark) is then 

(z 1 ) = (z 2 ) = l-±(z) c . 
It is convenient to write (z)c in the form 

(z) c = C F ^ {z + Cc l Z c } , 



where the universal C-independent piece is 



Z 



2 r 



V JO 



dz (3{z) [(1 - z - 2 7 )£(z) - (1 - z)\ 



The "hard" term Zq in Eq. (|45|) is given by 

Zy 



z, 



Ic-h, 

3-v 2 



with 



Ic 
h 



1 

V Jo 

1 

V Jo 



dz z 2 P(z) C{z) 



dz z 2 (3{z) . 



V Jo 



dz z 2 In 



l-(3(z) ' 



Performing the z integration in Eqs. ( p6j , |38D , we obtain 



Z 



U(3 + 2v 2 + 3v 4 )£ v -6(l + v 2 )} , 



(49) 




(50) 



(51) 



with 



C v = — In 

v 



l + v 



(52) 



1 -v 



The average gluon energy fractions for vector and axial production currents are 
obtained by combining the results of Eqs. ( f4l^|5"2"|) . Figure 2 shows the resulting (z) v 
and (z)a as a function of the top quark velocity v, for a s = 0.1. Note that in the 
relativistic limit, v — *■ 1, the average gluon energies diverge logarithmically, 



reflecting the emergence of the collinear singularity in the massless quark limit. Of 
course when a s ln(l/(l — v)) = 0(1) multiple gluon emission becomes important, and 
the leading logarithms must be resummed to all orders (see Ref. ||). The behaviour 
of the average energy in the threshold region, v — > 0, is discussed in the following 
section. 

Figure 3 shows that average gluon energy (E g ) = (z)\fs/2 as a function of the top 
mass M in e + e~ annihilation at centre-of-mass energy y/s. The correct proportions of 
vector and axial current contributions (Fig. 1) are included, and the strong coupling is 
taken to be a s = 0.1. For top masses in the experimentally favoured range 120 — 180 
GeV, (Eg) varies from 14 GeV to 4 GeV. 

2.5 Non-relativistic case 

For v 2 < 1, we can take the small- v limit of the results of the previous section to 
obtain 



(z) v w (z) A ^C F — \\ In — ^ 

7T Lo 1 — 



221 



9 J 



(53) 



Z 



16 



v 4 l + 0(v 2 ) 



(54) 



15 




(55) 



(56) 
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The high powers of the threshold factor v originate as follows: in the case of the 
"classical" universal term Z, one power of v 2 comes from the size of the integration 
region, z max = v 2 (Eq. (^)), and another power comes from the dipole suppression of 
the accompanying radiation, see also Refs. []T5|, [21|j . Because of the explicit z 2 factor 
appearing in the "hard radiation" spectrum, Eq. (^), the non-universal contribution 
Z c ~ 2^ ax = v 6 , and thus acquires an additional v 2 suppression relative to the univer- 
sal Z piece. Note, however, that the impact of the "hard radiation" in the threshold 
region strongly depends on the production channel. In particular, the Born factor 
in the denominator in Eq. (|45| ) can compensate the suppression of the non-classical 
effects in Zc, see Eq. ( |I2| , |T3|) . This is an interesting example of how the universal na- 
ture of soft radiation can sometimes be obscured by process-dependent short-distance 
effects. 

In the vector channel, where the (S'-wave) Born amplitude is not suppressed at 
threshold, the soft piece dominates and we find (see also fl5| , |2T|) 



(^>y~CV-^ 4 \l + 0(v 2 )} . (57) 

7T it) L J 

However, in the case of the axial- vector (P-wave) production channel, the non-radiative 
amplitude for the C — > qq decay is additionally suppressed at threshold. As discussed 
above, "non-classical" radiation effects at 0(a s ) are then equally important, and in- 
duce a significant contribution to (z)a- From Eqs. P5]j^j5^ - |5BD we find 

(z) A ~C F ^^v* (l + l) \l + 0(v 2 )} , (58) 



tt 15 V 7 

with the second term in the () brackets coming from the "hard radiation" . The limiting 
v — > behaviours corresponding to the leading v 4 terms in (z)y and (z)a are shown 
as dashed lines in Fig. 2. 

As we have already explained, the practical consequences of the differences be- 
tween the average gluon energy in the vector and axial-vector channels for e + e~ — > tt 
are not large, since the axial contribution is so small in the threshold region, see 
Eqs. ( ITT , VI , rBLIITI ) . It is, however, worth mentioning that the non-universality of the 



soft radiation in the threshold region is not restricted to the axial channel. The same 
situation arises in the case of tt production in the scalar (e.g. Higgs exchange) chan- 
nel, S — > qq. In the Appendix we present the corresponding results for (z)s- We find 
that in the limit v — > we have (z)$ — (z)a, so that once again the soft radiation 
contribution is accompanied by a hard radiation term of the same order. 

Another example which could be of practical interest concerns gluon radiation 
accompanying two scalar quark (e.g. stop) production: 

e + e~ — ► qq . (59) 
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This is a P-wave process, dominated by photon exchange.^ Near threshold, the average 
energy loss is again given by Eq. (|58|), illustrating the non-universality of the soft 
radiation description for this process also. 

All the results derived above have assumed that the top is a stable particle. However 
in practice and particularly near threshold, one must be careful to take account of the 
effects of the top decay width and of non-perturbative fragmentation. Although the 
former can affect the distribution of soft gluons (i.e. gluons with energy E g < Y t ) 
0, [| , it is not important for the average top energy loss. We can quantify the possible 
impact of non-perturbative dynamics by using the string-model results of Ref. |29| 



for the average amount of energy lost by the top quark before it decays. Near the tt 
threshold this leads to 

with the so-called string tension factor k ~ 1 GeV/fm. Since T t ~ M 3 ||, these 
non-perturbative effects decrease as M~ 4 . Taking as a specific example the canonical 
values M = 150 GeV and r t = 0.8 GeV, we obtain 

(A^) non pt ~ 3.3 x 10" 3 v . (61) 

Comparing Eqs. (|57| ) and (|6~l]), one may conclude that the perturbative prediction is 
only valid for v 2 > 0.17. 



3 Conclusions 

It is important that gluon emission in top quark production processes is under control. 
Otherwise, it may provide a serious source of uncertainty in the determination of the 
top quark parameters from experimental data. In this paper we have studied gluon 
bremsstrahlung off the top quarks in e + e~ — ► tt. We have calculated the exact first- 
order perturbative expressions for the energy fraction distributions of the i-quarks. 
These are valid in a wide energy range above the tt threshold. 

To examine the impact of gluon emission on top quark distributions we have cal- 
culated the average energy loss at the production stage. This infra-red safe quantity 
characterizes, for example, the difference between the beam energy and the actual 
energy shared by the produced t-quarks. 
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12 



We are grateful to Victor Fadin, David Miller, Lynne Orr and Torbjorn Sjostrand for 
useful discussions. This work was supported in part by the United Kingdom Science 
and Engineering Research Council. 



13 



Appendix 

In this Appendix we present the expressions corresponding to single primary gluon 
emission in the process e + e~ — ► qq, mediated by the exchange of scalar S or pseu- 
doscalar P particles (for example, the Higgs bosons of the Standard Model or its su- 
persymmetric extensions). Now the invariant polarization operator bears no Lorentz 
indices, and we define (cf. Eq. (|25|)) 



n c = I E Tr f(Pi + M)A C (M - p 2 )A c ] , (Al) 



A 



where the normalization has been fixed to preserve the relation fl30|). Here the scalar 
functions A c are analogous to of Eq. fl2"0|) but with the vertex operators 



75 , (A2) 



substituted for T^. 



Defining the double-differential distribution analogously to Eqs. (p£|), (p0|) with 

Cs = (a = v 2 , (p = 1 , (A3) 



we find 



n * - 4'" 2 + KS+S+ 2 i- (A4) 

. 2 q 2 1 f K 2 , Kl 



Hp = - - + - + 2 . (A5) 

The gluon energy spectrum is then given by Eq. ( p9f) with 

£h(z,S)=Z h (z;P) = £(z) , (A6) 

and the mean gluon energy fraction by Eq. ([45]) with 

Z s = Z P = I c . (A7) 

Near the threshold for the pseudoscalar channel, the soft term dominates and 

(z) P ~ (z) v . (A8) 

In contrast, the transition + — > qq is P-wave, and the Born term (s — v 2 compensates 
the additional suppression of Zs- The threshold result for (z)s therefore coincides with 
that for (z)a, Eqs. 
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Figure Captions 



[1] The functions t v (s) and r (s) denned in Eqs. (14,15) 



[2] The average gluon energy in e + e~ —>■ tig as a function of the top quark velocity 
v, for vector and axial vector current production. Also shown (dashed lines) are 
the limiting v — > behaviours given in Eqs. (|57j]58D , and the limiting v — > 1 
behaviour given in Eq. 

[3] The average gluon energy in tt production in e + e~ annihilation at t/s = 500 GeV, 
as a function of the top quark mass M. The strong coupling is fxed at a s = 0.1. 
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Fig. 1 
Fig. 2 
Fig. 3 
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